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Abstract 

We perform Monte Carlo simulations of a lattice gauge system with an action which 
contains two parallel plaquettes. The action is defined as a product of gauge group 
variables over two parallel plaquettes belonging to a given three-dimensional cube. The 
peculiar property of this system is that it has strong degeneracy of the vacuum state 
inherited from corresponding gonihedric Z2 gauge spin system. These vacuua are well 
separated and can not be connected by a gauge transformation. We measure different 
observables in these vacuua and compare their properties. 
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1 Introduction 



In this article we will report some observation that we have made in gauge invariant model 
defined on a lattice with an action which contains two parallel plaquettes . This action is 
non-Abelian generalization of the Z2 gauge spin system which was introduced as a lattice 
realization of random surfaces with gonihedric action It is defined as a product of 
gauge group variables over two parallel plaquettes belonging to a given three-dimensional 
cube of the lattice 

Hgonihedric = — ~ i^/n'^)Re TvUplaq " Re Trllplaq}, (1) 

^ {parallel plaquettes} 

where Upiaq is a product of gauge group elements Uij defined on a link < i,j > over all 
sites of a plaquette and Uij = exp{igA^a). As one can see the action contains a product 
of two traces taken over two parallel plaquettes belonging to a 3D-cube. The summation 
is over all pairs of parallel plaquettes 0. 

It is easy to see that for smooth classical fields the action reduces to the original action 
for non-Abelian fields as it takes place for the Wilson action 0, 

2r? 

H=— {l-{lln)ReTrUpiaq). (2) 

^ {plaquettes} 

Indeed the matrix product of four U' s around plaquette = 1,2 becomes 

Re Trllplaq = Re Tr exp{iga?Fi2 + ...) 
and after expanding fields around the center of the 3D-cube we will get 

1 - {l/n'){n - '^TrFUx, + ^S,Mn - '^TrFUx, + ^5,s)) - '-^TrF?2 

therefore it reproduces the classical action for pure Yang-Mills fields. 

Thus on a classical level both models (|l|) and (H) are equivalent. But on a quantum 
level it is not obvious at all. The reason is that the vacuum structure of the model (0) 
is drastically different 0, [^l- In addition to the trivial vacuum configuration when all 
plaquette variables are equal to one 

^Tr Uplaq = +1, (3) 

there are 2'^^ different vacuum configurations with frustrated plaquettes |Tr Upiaq = 
— 1. One of these vacuua can be easily described as a configuration with all plaquettes 
frustrated 

^Tr Uplaq = -1. (4) 

It is obvious that it is impossible to connect vacuum configurations (H) and (Q) by a gauge 
transformation. The riach vacuum structure inherited from the "parent" Z2 gauge spin 
system with the Hamiltonian [Q] 

Hgonihedric = ^ Y {(Tcxaa) {(Tcrcra) (5) 

^ {parallel plaquettes} 



■^The model with an action which is defined as a sum of contributions over all closed loops made up 
with six non-repeated links has been considered in Iq]. 
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In both models (|^) and (|1]) we have the same vacuum structure: towers of frustrated 
plaquettes describe different vacuua 0]. The difference between SU{2) and Z2 systems 
is that in the model (|ip we have a continuous Lie gauge group and therefore continuous 
fluctuations of gauge fields around these vacuua. 

The problem of our main concern is the physics in these vacuua and its relevance to 
the continuum limit. It is also interesting to understand: what is the remnant of this riach 
vacuum structure in the continuum limit. Is there any dilaton-like field which describes 
different vacuua in the continuum limit or not? And if there is a continuum limit in 
all these vacuua, then what is the difference between them? In the next section we will 
present results of simple analytical consideration of the model and the results of the Monte 
Carlo simulations for the SU{2) gauge group which are performed on three-dimensional 
lattices for an obvious simplification 0. 



2 High and low temperature expansion of average 
action 

As usually, all observables can be constructed as a product of gauge variables over closed 
loops and we shall consider first the average action. For the Wilson action (|^) a simple 
one-plaquette average is equal to [0 

P = < 1 - -Re Tr Upiag > = 1 - ^ + ... when p^Q 

and 

3 

P = — + ... when 3 CO. 

Let us now consider two-plaquette action (|l]) which, as we shall see, has a different 
high temperature expansion. At high temperatures the average action is equal to: 

PP = < 1 - Tr Upiag ■ Re Tr Upia, > = 1 - ^ + ... 
and the low temperature expansion is the same as for the standard action 

On Figure 1 one can see the behavior of average action in both models which we got by 



Monte Carlo simulations. In the simulations we used the Heat Bath algorithm |TT|. The 
new action has indeed smaller values at high temperatures and the slope is four times 
smaller. 

One can see also crossover from high to low temperatures. In the thermal cycles we 
have seen small hysteresis loop, but the measurement of susceptibility demonstrates that 
the maximum does not scale with the volume ( see Figure 2,3). Thus we do not see any 
phase transition at /5 ~ 3.5. 
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SU(2) phase diagram in 3-d for two different actions 
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Figure 1: The behaviour of the average action in two different models. The new action 
has four times smaller values at high temperatures. 
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Figure 2: The thermal cycle of the two-plaquette action. 
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Figure 3: Susceptibility versus (3 for two volumes 8^,24^. The maximum of susceptibility 
does not scale with the volume. 



3 Loop averages and vacuum structure 

In the model with a one-plaquette action (^) the high temperature expansion of the loop 
averages is equal to 



< W{C) > = W{I, J) 



ij 



and if one defines W{I, J) ~ exp{—aA), where A is the area IJ of the loop C then the 
string tension is equal to 



1 , A. 



In the two-plaquette model (|l|) the loop averages demonstrate peculiar behaviour. 
Indeed at high temperature the loop average 



< W{C) >=Z- ■ J dU e'^ ^ Tr. 



c Uij 



(6) 



is equal to zero for SU{2) gauge group, because there are no low order lattice diagrams 
(proportional, let us say, to area) which can contribute to the expectation value of W{C) 
(see Figures 4 and 5). We have absolute confinement in this model. It is similar to the 
limit 77, — > cxD of the model with one-plaquette action 



1 , 2n^ 



oo 



and thus < W{C) >=^ exp{-aA) 00. 



■^The situation changes for large gauge groups where there exist high temperature diagrams which 
contribute to one-plaquette average. Specifically for the SU{3) group the 3D-box diagram on a top of 
plaquette gives nonzero contribution. 
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Figure 4: Schematic multivacuum structure in the model with two-plaquette action. This 
figure shows the absolute confinement at high temperatures /3 < 3.5 and gauge nonequiv- 
alent vacuua at low temperatures. 



The Monte-Carlo simulations of SU{2) group demonstrate that this behaviour takes 
place until (3 ~ 3.5 — 4 as one can see on Figures 4 and 5. At these temperatures there 
appears proliferation of large lattice diagrams which cover almost the whole lattice and 
< W{C) > gets nonzero value. 

At low temperatures the nontrivial vacuum structure begins to play a central role. To 
see that let us consider the situation when two opposite plaquettes inside a 3D cube are 
frustrated. As it is easy to see from the action the energy of the new configuration will 
be the same. Indeed the energy is minimal when \TrUpiaq ■ TrUpf^^^ = 1, therefore we 
have two possibilities ether 

It^.. TT _ ^rj, jjparal _ 

or 

1 1 

^Tr Uplaq = -Tr U^f^g^ = -1. 

On the whole lattice one can build towers of frustrated plaquettes so that the new con- 
figuration will have the same energy as the trivial vacuum |Tr Upiaq = +1. Thus the 
vacuum has global degeneracy as it was in the corresponding Z2 gauge spin system 0. 
This actually means that we have many vacuua in which the gauge field is very far 
from perturbative value. 

This complicated vacuum structure is refiected on Figures 4 and 5 which show multi- 
vacuum structure at low temperature. The top curve corresponds to the vacuum in which 
all plaquettes are equal to plus one, |Tr Upiaq = 1, and the bottom curve corresponds to 
the opposite case when all plaquettes are frustrated |Tr Upiaq = — 1- It is also not difficult 
to construct explicitly vacuum configurations in which only part of the lattice plaquettes 
are frustrated. The vacuum with all plaquettes frustrated, |Tr Upiaq = — 1, corresponds 
to a gauge field which is far from the perturbative value = in the whole lattice. 
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mean value for Wilson loops in the two plaquettes action, volume=8'*3 
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Figure 5: The behaviour of the mean value for the 1x1 and 2x2 loops for two-plaquette 
action. It represents the evolution of the averages starting from low temperature with 
|Tr Upiaq — +1 to high temperature and then back to low one. As one can see the system 
fluctuates at pseudocritical point and then falls into vacuum with Tr Upiaq ~ 



Two natural questions arize: 

i) is there a sensible continuum limit in all these vacuua in the limit a ^ 0, j3 ^ oo 7 

ii) and if there is a continuum limit, then what is difference between them ? What type 
of parameter ("couphng constant") or field ("dilaton") may describe them? 



4 String tension in different vacuua 

In this paragraph wc will study the behaviour of string tension in the limit (3 ^ oo. We 
choose this quantity because it can characterize the system in the continuum limit if it 
exists. As usually we define the string tension by the behaviour of W(L,L) for large L as: 

W{L L) Q-i'^L'^+mL+Const) 

To calculate string tension we shall follow the procedure: 

i) calculate W(L,L) for different values of L < 6 at a given /3 (in our simulations the 
lattice volume is 16^ and the coupling constant is in the range 6 < /3 < 10), 

ii) then fit the quantity —ln\W{L,L)\ using second order polynomials and extract a for 

each f3, 

iii) because string tension scales as a = aphygO? and (3 = 2/g'^a we shall plot the ratio 
Py/a^ versus 1//? (see Figure 6 ). 

The last ratio converges to a finite limit when (3 oo. We expect that this limit can 
be reached within the finite size corrections which arc of order O(^). On Figure 6 we also 
depicted string tension in the vacuum < |Tr Upiaq >= +1 calculated for the larger lattice 
24^ in order to see the variation with volume. The convergence is better for small (3. 
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Figure 6: String tension versus inverse [3 in the standard vacuum Tr Upiaq = +1. 

Until now we have been calculating string tension in the vacuum < |Tr Upiaq >= 
+ 1 when /? — > oo. In the new vacuua we have frustrated plaquettes and the quantity 
|Tr Upiaq is less than one. It remains constant during Monte Carlo simulations for large 
P >6 and we conclude that these vacuua are well separated 0. 

To explore vacuum structure we have to calculate W{L, L) for these vacuua as well. 
A typical result is shown in Table 1, where we choose /3 = 10, the volume is equal to 16^ 
and L is smaller than six. By "vacuum —1/3" we mean the vacuum which has the mean 
value < \Tr Upiaq >= -1/3. 

We observe that the absolute value |iy(L, L)| is the same for two vacuua +1 and -1, 
as well for the vacuua +1/3 and -1/3. The absolute value L)| for odd L in the 

vacuua +1/3 or —1/3 is one third of its value in the vacuua +1 or —1. Because of this 
relation, \Wx[L,L)\ = HWi/aiL, L)\ we have Q 

-ln\Wi/3iL, L)\ = -ln\Wi{L, L)\ + ln{3) = al'^ + ml + ln{3) + Const 

and the string tension is the same for all these vacuua. Thus we can argue that string 
tension does indeed scale in all these vacuua and we can speculate that one can distinguish 
them by the scalar operator | < Tr Upiaq >-gluon like condensate |12|. 

One of us (K.F.) wishes to thank the EU for partial financial support (TMR project 
FMRX-CT97-0122). 

TABLE 1 

The values of < W{L, L) > in different vacuua for f3 = 10 and volume 16'^ : 



^When we start from a given vacuum configuration < ^Tr Upiaq >7^ +1 it never drifts into other 
vacuua when /? > 6. 
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By W-i{L, L) we mean the mean value of W(L,L) in the vacuum Tr Upiaq — — ^ for example. 
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L 


vacuum +1 


vacuum -1 


vacuum -1-1/3 


vacuum- 1/3 


1 


0.946281(4) 


-0.946277(3) 


0.315429(3) 


-0.315427(4) 


2 


0.83866(2) 


0.83863(2) 


0.83865(1) 


0.83864(2) 


3 


0.71894(4) 


-0.71889(4) 


0.23968(3) 


-0.23967(4) 


4 


0.60273(7) 


0.60265(8) 


0.60269(8) 


0.60265(9) 


5 


0.49675(13) 


-0.49656(12) 


0.16561(9) 


-0.16567(10) 


6 


0.40381(19) 


0.40352(20) 


0.40361(19) 


0.40369(21) 
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